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Abstract 

Usually, the relaxation times of a gas are estimated in the frame of 
the Boltzmann equation, which is a reduced description if compared 
with the microscopic Hamiltonian model, and, moreover, being a dis- 
sipative equation, already has an irreversible character. In this paper, 
instead, we deal with the relaxation problem in the frame of the dy- 
namical theory of Hamiltonian systems, in which the definition itself 
of a relaxation time is an open question. We introduce a suitable def- 
inition of relaxation time, and give a general theorem for estimating 
it. Furthermore we give an application in a concrete model of a gas, 
providing a lower bound to the relaxation time, which turns out to be 
compatible with those observed in dilute gases. 



1 Introduction 

The definition and the estimate of relaxation times are problems of central 
interest when one attempts at describing macroscopic systems through mi- 
croscopic Hamiltonian models. 

In the case of gases, these problem are tackled, and solved, in the frame 
of the Boltzmann equation (see [1]). In such a frame the existence of a 
relaxation time is somehow obvious, due to the irreversible character of the 
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equation, and the estimate is obtained in terms of the eigenvalues of the 
hnearized equation, about the equihbrium solution. 

There are two problems left. First, one doesn't know how to proceed in 
the case of a system different from a gas. Furthermore, even in the case of a 
gas, the Boltzmann equation refers to a reduced description, whose relation 
to the original Hamiltonian system is not completely clarified. 

An approach to the relaxation problem that refers to the complete system, 
without any reduction, was given recently (see [2]). A characteristic feature 
of such an approach is that the system is assumed to be non isolated, being 
in contact with a mechanical thermostat which makes the system dissipative 
albeit reversible. 

In the present work we tackle the problem from the point of view of the 
dynamical theory of Hamiltonian systems, for systems which are isolated. In 
this perspective, a partial answer to the problem is given by Kubo's linear 
response theory [3]. Indeed, such a theory enables one, at least in principle, 
to compute in microscopic terms the macroscopic transport coefficients, and 
then, via macroscopic equations, the relaxation time. From our point of 
view, however, this answer is not completely satisfactory, because it appeals 
to macroscopic irreversible equations, which should preliminarly be deduced 
from the microscopic ones. 

A related but different approach is followed here, whose main scheme can 
be sketchted as follows. 

From linear response theory we take the starting point, namely, the idea of 
following the time evolution of the probability distribution function in phase 
space (and not in the reduced /^-space, as in the Boltzmann equation), when 
a perturbation —hA{p, q) to the original Hamiltonian Ho{p, q) is introduced. 
Still following Kubo, we then choose to concentrate our attention on a par- 
ticular observable, namely the one conjugated to the perturbing field, in the 
familiar sense in which pressure is conjugated to volume and magnetization 
to the magnetic field. As an example, later on in this paper we will deal with 
the simple case in which the perturbing field is gravity, and the conjugate 
observable is the height of the center of mass. However, while Kubo looks 
(in a heuristic way) at the asymptotic behaviour of the distribution function 
(and hence also at the expextation of the conjugate variable) for large times, 
when equilibrium has been attained, our attention is instead addressed at 
defining and estimating the relaxation time itself. 

In the spirit of the Kubo approach it is natural to say that equilibrium 
is attained when the time derivative of the conjugate variable is negligible. 
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The main idea is then that the relaxation time should be defined as the one 
at which the rate of change of the conjugate variable is smaller than a certain 
threshold. 

The aim of the present paper is indeed to give a lower bound to the 
relaxation time, looking at the evolution of the time-derivative of the variable 
conjugated to the perturbation A in the Hamiltonian. It is easily seen that 
the time-derivative of the conjugate variable is the function B =^ [Hq, A], so 
that this is the quantity on which we will concentrate in this paper. Having 
chosen the relevant function, namely B, we make use of the easily established 
properties that its expectation vanishes at equilibrium, and that its time- 
derivative is positive at the initial time. Thus a lower bound treiax to the 
relaxation time is provided by the time before which the time-derivative of 
the expectation of B is proven to be positive. 

The problem is then that one should make use of suitable a priori esti- 
mates on the dynamics, in lack of an explicit integration of the equations 
of motion. This can actually be implemented following the main idea in- 
troduced in paper [1], which was concerned with Hamiltonian perturbation 
theory in the thermodynamic limit. In such a paper, a procedure is given 
which, for any function / of phase space with respect to Gibbs measure, 
allows one to provide an upper bound to \\Utf — /II2, by knowing an upper 
bound of ||[/, if]||2. Here, [■, ■] denotes Poisson bracket, H the Hamiltonian 
of the system, and Ut the corresponing unitary evolution group. 

The estimate of the lower bound tj-eiax is provided by formula (Q of The- 
orem [H which is stated and proved in Section 2. Such a proof is given for an 
ample class of Hamiltonian systems, which are the ones considered in most 
rigorous works in Statistical Mechanics (see [5J). 

In Section 3 the general theorem is applied to the case of a gas of non- 
interacting point-particles enclosed in a cubic box, to which the gravity force 
is added as a perturbation. The relaxation time thus found turns out to be 
comparable with the typical relaxation times observed in dilute gases. One 
should actually study the model of a gas of particles interacting through a 
stable and temperate two-body potential. Working at a heuristic level, it is 
easily seen that, for a dilute gas, essentially the same estimates are obtained 
than in the case of a non-interacting gas. We hope to be able to give a 
rigorous treatment in the future. 

Some further comments are given in Section 4. 
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2 General theorem about relaxation times 



Consider a system, with phase space M., which is in equihbrium with respect 
to the Gibbs measure relative to an unperturbed Hamiltonian Hq at inverse 
temperature /3, i.e., the measure with density po given by 

Po = ^ exp {-(3Hq) , 

Zq being the partition function. Suppose at time a perturbation —hA{p, q) 
is introduced, where A is a given function on phase space and the parameter 
/i > controls the size of the perturbation. So, at positive times the Hamil- 
tonian is Hi = Ho — hA. The corresponding Gibbs density will be denoted by 
pi. Our aim is to find a sensible definition and an estimate for the relaxation 
time treiax to thc final equilibrium with respect to Hamiltonian Hi. 

To this end, along the scheme sketched in the introduction, following 
Kubo we consider the observable B defined by 

B=[A,Ho] , 

i.e., the time derivative of the perturbation A with respect to the flow gener- 
ated by the total Hamiltonian Hi. We then consider the probability density 
p, solution of Liouville's equation relative to the total Hamiltonian Hi with 
initial condition p{0) = pq, and look at the evolution of the expectation of 
B, i.e., we look at the quantity 

B{t) = J Bp{t)dpdq . 

The quantity of interest actually will be its increment 

AB{t) = B{t) - B{0) . 



Writing p in the form 



Pit) Po + Ap{t) , (1) 



by Liouville's equations the perturbation Ap satisfies the differential equation 
dAp 



dt 



[Ho-hA,Ap]-h[A,po] , (2) 
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with Ap(0) = as initial condition, and one has 



AB{t) = j B Ap{t)dpdq . (3) 

We will show that under the familiar conditions which entail reversibility 
(namely, that both Hq and A are even in the momenta), the quantity AB 
vanishes not only (as it is obvious) at time zero, but also at equilibrium 
with respect to the total Hamiltonian Hi. This is due to the fact that the 
expectations of B with respect to the Gibbs densities po and pi corresponding 
to the Hamiltonians Ho and Hi, both vanish by symmetry, because po and 
pi are even in the momenta, whereas B is odd. On the other hand, it 
turns out that AB is initially an increasing function of time, since its time- 
derivative is positive at time 0, as it will be shown later. Thus, the time- 
derivative of AB has to become negative if equilibrium with respect to the 
total Hamiltonian has to be attained, and consequently a lower bound to the 
relaxation time is provided by the time treiax up to which the time derivative 
of AB is guaranteed to be positive. 

We thus define the relaxation time treiax by 



treiax = SUp t* , (4) 



where t* is such that 



^A5(t) > for all < t < t* . (5) 

The problem is then to estimate the rate of growth of AB. Now, on the 
one hand, following Kubo we know that AB(t) is strictly related to the time- 
autocorrelation of B (see (1131) below). On the other hand, we can make use of 
the main result obtained in paper [1] , in which it was shown how to estimate 
the time-autocorrelation of B in terms of the Hamiltonian. Indeed, from the 
main result of such a paper one easily obtains the following property: an a 
priori estimate of the type 



\[B,Ho]\\,<v\\B\ 



' 



(with the norm defined below) implies that the time-evolution of AB is slow 
if rj is small, or, more precisely, that the relaxation time defined by (jlj), ([5]) 
is inversely proportional to rj. 
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Here, || ■ ||o is the norm on C^Ai), the Hilbert space of square integrable 
complex functions on A4, with respect to po- We will also have to consider 
the Hilbert space Cl{Ai) of the square integrable complex functions with 
respect pi. The corresponding £^-norm will be denoted by || ■ ||i. 

Under the rather natural condition ([6]) given below, which ensures the 
smallness of the "change" of the Gibbs measure induced by the perturbation 
—hA, for a large class of observables it can be proven that the two norms 
just introduced are asymptotically equivalent as h ^ 0. Indeed one has the 
following lemma, whose proof is deferred to Appendix A. 

Lemma 1 Assume there exist 5 > and K > such that 

[ dpdqe^^po<K and [ dp dq e-^^ po < K . (6) 
Jm Jm 

Then, for all real functions f on M. satisfying at least one of the conditions 
||/2||^<+oo, ||f||^<+oo, (7) 



one has 



^ " "^"g = o(l) , as h Q 



We are now able to give an estimate for t^eiax in terms of r], which is 
provided by the following Theorem 1. It will be seen that some technical 
hypotheses, namely those given in ([8]) below, are required just in order that 
at least one of the conditions ([7j) of Lemma [1] is satisfied. 

Theorem 1 Let the unperturbed Hamiltonian Ho{p,q) be even in the mo- 
menta and bounded from below, and consider a perturbation —hA{p,q), with 
h > and A even in the momenta. Suppose A and Hq are such that hypoth- 
esis ^ of Lemma 1 is satified. With B = [A,Hq], suppose furthermore that 
the following technical conditions are satisfied: 

\\B\<+cx; \\[B,Hof\l<+^ ; \\[B, A]\ < +cx^ . (8) 
Then a lower bound to the relaxation time defined by ^ is given by 

trelax > — + aS k ^ , (9) 

where rj is such that 

||[i?,//o]|lo<^l|i?llo • (10) 
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Remark. It may appear that some conditions are too restrictive for the 
Theorem to be used in the thermodynamic limit, but it turns out that such 
a difficulty can be overcome. For example, Hq was required to have a finite 
lower bound, call it D; however, the result is found not to depend on the value 
of D. So, D can grow with the number N of degrees of freedom, without 
affecting the validity of the Theorem, provided D is finite for any finite A^. 
A similar way of proceeding works also as far as conditions (jH]) and ([H]) are 
concerned, so the Theorem holds for any system, however large it may be. 
This fact allows one to pass to the thermodynamic limit. 

Proof of Theorem [T], First of all, we notice that equation ([2]), which 
governs the time evolution of the pertubation Ap, admits a unique solution 
in the Hilbert space Cl{Ai). Indeed, equation ([2]) is a linear inhomogeneous 
first-order differential equation in CKAi) of the form 

x = dx + f , 

where the operator O =^ [Hi, ■] generates a semigroup of unitary evolution 
transformations (see for example [H]). Thus, since the second term /z[A,po] 
at the r.h.s. belongs to £l{Ai), as will be shown below, the solution is known 
to exist and to be unique (see Theorem 3.3, page 104, of [7]). Such a solution 
is given by a simple adaptation of the variation of constants formula, namely 

Ap(x,t) =(3h [ ds5($'x)po($'x) , (11) 
Jo 

where x =^ {p,q) denotes a point of phase space Ai, and the fiow gen- 
erated by Hi. Notice that, as the initial datum vanishes, one obviously has 
OAp(O) e Cl{M). 

We show now that [A, po] G Cl{M.), too. To this end, we first notice that 

\\[APo]\\i = \\Bpo\\i<—\\B\\, , (12) 

^0 

def 

where D = inip ^Ho. On the other hand, iterating the Schwarz inequality 
gives 

il^ilo<(||5llo)^<(ll^1lo)'<+^' 
in which the first hypothesis of ([8]) was usecl^- In virtue of such an hypothesis, 

^According to the same reasoning, the square of the norm of a function will be bounded 
from above by the norm of the squared function. 
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we can also apply Lemma[T]to B and observe that, for h small enough, \\B\\-^ 
is finite. Thus, by (IT^ it is proved that [A,Po] belongs to CKAi). 

We now look at the expextation B{t) and at its increment AB{t). By 
using ffTTj) for Ap in ([3]), one finds for AB{t) the expression 



M 



AB{t) = l3h dx / c/s5(<l>^x)B(x)po($'x) 



(13) 



Using the shorthand /(x^) = /($ *x), one has then: 



^A5(t) = / dx5(x_i) fi(x) po(x-i; , 

M 



or equivalently (due to preservation of Lebesgue measure), 
^AB{t)=(5h I dx5(xi)5(x)po(x) . 



(14) 



M 



At this point we remark that the integral in f[T^ could be evaluated in 
a quite simple way, if there appeared pi in place of po- Indeed, due to the 
unitarity of the fiow, for any / in C\{M.) one would have 



dx/(x,)/(x)pi(x) 



M 



?-^ii/(x.)-/(x)ii; 



(15) 



and thus, on account of hypothesis ffTOl) of the Theorem, the thesis would 
follow by using Theorem 1 of [Ij (see below). 

The rest of the proof is devoted to show that the error made by taking 
Po in place of pi is negligible in the limit /i ^ 0. To this end, we suitably 
rewrite f[T^ in the form 



dt 



AB{t) =ph 



dx5(xt)5(x)pi(x) - / dxfi(xOfi(x) (pi(x) -po(x)) 



M 



>/3h 



dxS(x05(x)pi(x) 



M 



dxS(xOS(x) (pi(x)-po(x)) 



M 



First, we show that the second term at the r.h.s. vanishes as h 
Indeed, by Schwarz's inequality we have 



.(16) 
0. 



dx5(xi)5(x) (pi(x)-po(x)) 



M 



< 



dx (52(xi)52(x))pi(x) 



M 



(17) 
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where we have defined 



M 



This function coincides with the one defined by 0341) in Appendix A. As there 
shown, 7(/i) — > as /i — 0. We then make use of f|T5l) . by replacing for 
/ and neglecting the negative term, to find an upper bound to the r.h.s. of 
ffTTl): one has, in fact. 



dx (i?2(x05'(x))pi(x) 



M 



In order to show that ||-B^||i is finite, we use Lemma [H whose hypotheses are 
satisfied owing to the first inequality of ([HD- Thus, as ^{h) ^ for /i — 0, 
one gets 



dx5(x05(x) (pi(x)-po(x)) 



M 



o(l) as /i — > 



(18) 



We then come to the first term at the r.h.s. of (fT6l) . which, using (fT5l) 
again, can be estimated as 

/ dxi?(x,)fi(x)pi(x) = ||i?||?-i||5(x^)-i?(x)||; . (19) 

We now make use of Theorem 1 of [3], which ensures that, if 

\\{B,E^\\\^<^\\B\\^ 

is satisfied, then one has 

||5(xO-5(x)||, <r/t||5||, . (20) 

Now, to give an estimate for r/, we notice that the following inequalities hold 
as /i — 0: 

< ||[i?,/7o]|li + /^||[5,^]|li 

< (||[5,i7o]|lo + o(l))^ 
+mi[i?,A]||^ + o(l))^ . 
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Here, in the first line the triangle inequality was used, while the second line 
is a consequence of Lemma [H the hypotheses of which are satisfied in virtue 
of the second and the third inequalities in ([8]). Hence, by hypothesis ffTOj) we 
obtain 

\\[B,H,]\\, < \\[B,Ho]\\, + oil)\\B\\, 

< v\\B\\, + oil)\\B\\^ 

< {v + oil))\\B\\, , (21) 

so that fj = r] + o(l). 

Eventually, by replacing in (fT6l) the estimates (fT8|) and (fT9|) . one has 

AA5(t) > (^1 - !^ - o(l) - o(l) ■ t'^ \\B\\l . (22) 

Therefore, in the limit as h ^ 0, the time derivative of AB{t) remains 
positive for 

t < — + 0(1) 
V 

and this completes the proof. 

Q.E.D. 



3 A gas in a gravitational field 

We study now a gas of non-interacting particles enclosed in a tridimen- 
sional box of side L. A comment on the possible extension to the case of 
interacting particles will be made later. Our aim is to show that Theorem 1 
holds if we take as a simple example of a perturbation the force of gravity, 
in which case the conjugated variable will be the vertical component of the 
total momentum of the system. 

For what concerns the interaction of the particles with the walls, due to 
the form of the conjugated variable it turns out that only the interaction 
with the horizontal walls will matter. Thus we limit ourselves to choose a 
particular form for the interaction potential with the horizontal walls. The 
unperturbed Hamiltonian Hq is then 

N 2 N 

^o'=^Ey + E^(-.)' (23) 
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where G (—00, +00), G (— L/2, -L/2), a = 1, 2, 3, g| = and V denotes 
the interaction potential with the horizontal walls, which we take as the 
repulsive part of the Lennard- Jones potential, namely as 



V{z) = 6 



12 / \ 12 
1 \ 1 



(24) 



where 5 is a positive parameter. According to the general scheme previously 
discussed, we add at time t = a perturbation —hA, and for the observable 
A we make the choice 

N 

^ = E-- (25) 

We then have 

Theorem 2 Let Hq and A be given by [23^) and / f^) respectively. Let the 
parameters 5 and L in [2^) be such that {[35)^2 < where (3 is inverse 
temperature. Then, as h 0, for the relaxation time one has the estimate 

t relax > to, With 

to = -(/35)^V^, (26) 
c 

and 

J u^2e-^duj ^ 25 . (27) 

Remark 1. By inserting the proper dimension, condition (PS)^ < L/3 
turns into (/35)t2 a < L/3, where a is the characteristic parameter for the 
range of the Lennard- Jones potential. This just expresses the requirement 
that the interactions with the walls have a range which is negligible with 
respect to L. Notice that we inserted the factor 1/3 at the r.h.s. just in 
order to fix a numerical value for c, but any other reasonable choice would 
not affect the result. 

Remark 2. The time to, once the correct dimensional constants have 
been introduced, becomes 



tn 



-{P6)^^/PmLa 



where m is the mass (of molecular order) of each particle. Therefore, for 
macroscopic systems in which L is of the order of magnitude of 1 m, while 
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10~^°m, at room temperature, one gets to ~ 10 ^s, a value which is 
of the order of magnitude of the typical relaxation times measured in gases. 
Further comments will be given in the next Section. 

Proof. The proof consists in showing that the hypotheses of Theorem [T] 
hold, with 

V2c 
V = TTTTT 

' ^ '') 24 - 



In the first place, (Q is satisfied for any 6, because it involves integrals of 
continuous functions over a compact, and the integral over the p coordinates 
is equal to 1. 

As far as the other hypotheses are concerned, we first notice that in the 
present case one has 



The integral over the momenta of any power of B can be easily turned into 
a combination of terms of the form 



+ 00 



dpp'^e' 



which are finite for any n, thus proving the first of ([8]). In particular, one has 



\B\ 




(28) 



Clearly, one also has 

[B,A]=N , 

and hence the third of ([8]) holds. We then compute || [i?, iJo] llo- One has 



N 



13 



13 



Zj + 



+ 



(29) 



This function, and its powers too, are actually singular at some point in 
phase space, but they diverge there as a power, while the density po vanishes 
as an exponential, making the norm finite. So the second of (IHl) is proved, as 
well. 



12 



There finally remains the task of providing an estimate for the quantity 
rj, which is an upper bound for the ratio || [S, iJo||o/||-B||o- To this end, from 
(|29|) we get that \B, Hq]'^ is equal to 



14452^ 



13 



13 



13 



13 



^l + l 



and we have to estimate its po norm. Now, all terms with j 7^ / have vanishing 
integral, because they are products of independent quantities with vanishing 
expectations (being odd functions of z, while the density is even). So, there 
remain only the N terms in which j = I. Each of them contributes to the 
sum as 

26 



1445^ 
Z 

in which 



dz 



26 



1 



z + 



+ 



z — 



if 



-t35V{z) 



(30) 



dze-^'^^^^ . 



This integral cannot be computed in a closed form, so we give here an esti- 
mate: numerical computations show that Z is larger than L/4 if (/35)i2 < 
L/3. In the integral (!30|) . the third term is always negative, while the two 
other ones are identical, due to simmetry. An upper bound to each of them 
is given by 

\ 26 

1 



dz 



-l35V{z) 



< 



< 



dz- 



(35_ 
ZT2 



(/?5)i JO 



+00 



A quick comparison with 



then leads to 
V2c 



ojllo 



< 



\B\ 



' 



(31) 



with c defined in ( l27l) . 

This relation, on account of Theorem [H leads to formula ( 126 



Q.E.D. 
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4 Conclusions 



We have provided by Theorem [T] a lower bound to the relaxation times in 
Hamiltonian systems, and shown that in the case of a non-interacting gas 
enclosed in a box such an estimate is of the order of magnitude of the typical 
relaxation times measured in gases. 

Obviously, the interesting case to be studied would be that of a gas of 
interacting particles. At the moment, we are unable to provide rigorous esti- 
mates for such a case, so we decided to limit ourselves to the non-interacting 
case, and we leave for future work the interacting case. The following com- 
ment is however in order. If two-body internal forces were present, it is 
easily seen that their contribution to [B, Hq] (and to B itself) would vanish. 
The only difference would concern the estimate of the norm of the expression 
(!29|) for [B,Ho], but a heuristic argument indicates that, for sufficiently low 
densities, the difference with respect to the non-interacting case would be 
small. This argument becomes a rigorous theorem (see [8]) in the thermo- 
dynamic limit, for forces due to stable and temperate two-body potentials. 
So, in order to give such an estimate for finite N and compute the integrals, 
we would need at least an estimate of the convergence rate of the Mayer- 
MontroU equation, in the canonical ensemble, which is not known to us. On 
the other hand, this fact seems to indicate that the interactions with the 
walls, which we have considered in the present work, might have sensible 
effects even when one is interested in investigating relaxations of observables 
related to internal interactions. 

We now add some comments concerning possible further developments. 

The first point concerns the hypothesis made in Theorem [H that Hq and 
A are even in the momenta, which actually is not at all essential. Indeed, if 
such an hypothesis is not satisfied, it suffices to define Ueiax in a different way, 
namely, as the time up to which the time-derivative of AB remains larger 
than, for example, |^Ai?(0). This way, by inequality fl22l) one could prove 
Theorem [H except for setting l/?7 in place of V^/rj in ([9]). We decided to deal 
with the case of reversible Hamiltonians just because it is a very important 
one; furthermore, in such the relaxation time can be defined with no 

reference to arbitrary features, as the factor 1/2 introduced above. 

As a more interesting fact, we are confident that our line of reasoning 
may be extended to the case of perturbations of a finite size h, because this 
would just entail to consider the norms in Cl{M) rather than in Cl{M). 
Indeed, if we substitute [B, Hi] for [B, Hq] in hypothesis ffTOjl and use there 
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the norm || ■ ||i instead of || ■ ||o, we can directly set tj — tj in fl2UI) and there 
is no need of deducing (pTj) . so that the second and third conditions in ([H]) 
are no more required. Moreover, the first condition in ([8]) can be replaced by 
the condition that is finite, which makes trivial the proof that [A, Hq] 

is in Cl{M). 

The really open problem that remains in order to implement an extension 
to the case of finite h, at least for macroscopic systems, is the estimate of 
the difference of the two norms in (ITTj) . The estimate which appears in such 
a formula has the serious flaw of increasing exponentially with the number 
N of particles. This occurs because the upper bound provided there, which 
is an immediate consequence of Schwarz inequality, is valid for all functions 
in A way to improve such an estimate would be to restrict oneself to 
perturbations having some suitable characteristic feature. 

In particular, the work [9j of Lanford seems to suggest a good starting 
point. There it is pointed out that the only observables of interest in de- 
scribing a macroscopic system are the ones he calls finite range observable^, 
namely, observables which are sums of terms depending only on the position 
of a finite number of particles. The difference of the two norms in question 
can then be evaluated for each term and this should lead to an estimate 
which doesn't increase too much with N. We think that, if one limits one- 
self to considering a smaller class of functions, there is good chance that the 
problem of the number of degrees of freedom is overcome, and that some 
results are obtained also for the case of perturbations of finite size. These 
interesting investigations are left for possible future works. 

A Appendix 

Proof of Lemma [T]. We take as starting point the obvious equality 

/ (ipdqfpi= dpdqfpo+ dpdqf{pi-po) 
Jm Jm Jm 

which, by Schwarz inequality, gives 



/ dpdqf{pi-po) 
Jm 



< (^J^dpdqfpoyVW), (32) 



^As a matter of fact, he explains that this definition is chosen to make things simpler 
and is too restrictive. He gives also a reference to Ruelle's book in which it is shown 
how to deal with a broader class of observables, which represents the class of real interest. 
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with 



One can also write 



def 



dpdq 



M 



l{h) 



/ dpdge 



Pi 
Po 

2hl3A 



Po • 



Po 



. 1 , (33) 

as is seen by expanding the square and using the fact that po and pi are the 
densities of the Gibbs measures corresponding to Hq and ifi, respectively. It 
is also of interest to provide an upper bound to the l.h.s. of (132|) in terms of 
pi rather than of po- Indeed, one has 



dp dq f (p 



Po, 



M 



< 



dpdqfpi] \/^{h) 



M 



where, in a way similar to fl33|) . one gets 

def 



dpdq 



M 



Po 

Pi 

dp dg e^^^pQ 



2 

l] Pi 



M 



dpdqe-^^^po ) - 1 



M 



Now, we observe that the functions 7(/i) and 7(/i) can take arbitrarily 
small values as h goes to 0, if (E]) is satisfied. Indeed, by their definitions, 
they are always non-negative quantities. Thus, in order to show, for example, 
that 7(/i) < e for any fixed positive e, it will suffice that one has 

dpdq e^^l^^' 



Po 

(Xv(dpdge^^>o)^ 
To this end, let us note that 

1 



<l + e . 



/vjdpdge^^^Po 



< 



dp dq e 



-h/BA 



M 



according to Schwarz inequality. We combine this estimate with the Holder 
inequality, on whose account, ii h < ^, one has 



dpdge 



±chl3A 



Po 



< 



M 



dp dq e^^'^po 



M 

< K— 



chl3 
S 



dpdqpo 



M 



(34) 
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and we eventually obtain 



M / \JM 

< K— 1 as /i ^ . 



The immediate consequence is that, if we take h < min ^^^^^j, then 

'y{h) is less than e. An analogous argument, still based on inequality flM|) . 
ensures that 'y{h) takes arbitrarily small values as h goes to 0, too. Hence, 
the difference between and ||/||^ vanishes with h, provided p belongs 
to Cl{M) or Cl{M). 



Q.E.D 
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